The Pleiades mass function : models versus observations 
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Abstract. Two stellar-dynamical models of binary-rich embedded proto-Orion-Nebula-type clusters that evolve 
to Pleiades-like clusters are studied with an emphasis on comparing the stellar mass function with observational 
constraints. By the age of the Pleiades (about 100 Myr) both models show a similar degree of mass segregation 
which also agrees with observational constraints. This thus indicates that the Pleiades is well relaxed and that it 
is suffering from severe amnesia. It is found that the initial mass function (IMF) must have been indistinguishable 
from the standard or Galactic-field IMF for stars with mass m ;$ 2 Mq , provided the Pleiades precursor had a 
central density of about 10** stars/pc'^. A denser model with 10^'* stars/pc"^ also leads to reasonable agreement 
with observational constraints, but owing to the shorter relaxation time of the embedded cluster it evolves through 
energy equipartition to a mass-segregated condition just prior to residual-gas expulsion. This model consequently 
looses preferentially low- mass stars and brown dwarfs (BDs), but the effect is not very pronounced. The empirical 
data indicate that the Pleiades IMF may have been steeper than Salpeter for stars with m^2 Mq. 
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1. Introduction 

Because of its richness and proximity ('^ 120 pc, Robichon 
et al. 1999), the Pleiades cluster with an age of 120 Myr 
(Stauffer et al. 1998; Martm et al. 1998) is one of the 
best studied young open clusters. It has a mass of about 
740 M0, contains ~1000 stars and has a tidal radius of 
about 13.1 pc, a half mass radius of about 3.66 pc (Pinfield 
et al. 1998), and a binary-rich population with properties 
similar to that in the Galactic field (Bouvier et al. 1997). 
Recent stellar-dynamical computations suggest that the 
Pleiades-precursor may have been very similar to the 
about 1 Myr old Orion Nebula cluster (ONC) (Kroupa, 
Aarseth & Hurley 2001, hereinafter KAH), opening the 
possibility of seeing the same kind of object at two very 
different evolutionary stages. 

Several surveys have been performed to date to measure 
its stellar mass function which is now well constrained 
from the stellar to the sub-stellar domain. Recent esti- 
mates extend as deep as O.O3M0 (Moraux et al. 2003). A 
pressing issue in star formation is to know whether this 
mass function (MP) observed at an age of ^100 Myr is 
representative of the initial mass function (IMF). In other 
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words, is the observed Pleiades population representative 
of the cluster population at the time it formed and how 
did it evolve within 100 Myr ? 

KAH present calculations of the formation of Galactic 
clusters using Aarseth's Nbody6 variant GasEx (Aarseth 
1999). They defined an initial model from a set of assump- 
tions describing the outcome of the star formation pro- 
cess and studied how it could lead to the formation of the 
ONC and then to the Pleiades by dynamical and stellar 
evolution. They found encouraging consistency with ob- 
servational constraints on the radial density profile, kine- 
matics and binary properties for both the ONC and the 
Pleiades. With this contribution we compare the MF pre- 
dicted by the models with currently available data to test 
if the Pleiades IMF may have been compatible with the 
"standard" IMF (Kroupa 2001) or with the log-normal 
IMF found to match the Galactic field population down 
to the substellar regime (Chabrier 2002). 

In this paper, we describe briefly KAH's models in sec- 
tion 2 before comparing them to current estimates of the 
Pleiades MF and presenting relevant results about the ra- 
dial distribution and the evolution of the cluster popula- 
tion in section 3. The conclusions are given in section 4. 



I. Model description and predictions 

KAH develop stellar-dynamical models of cluster forma- 
tion under realistic conditions. Their computations are 
performed with Aarseth's NBODY6-variant GasEx which 
allows accurate treatment of close encounters without soft- 
ening and multiple stellar systems in clusters, incorpo- 
rates stellar evolution and a Galactic tidal field in the 
solar neighbourhood. In the models the initial stellar and 
sub-stellar masses are distributed following the standard, 
or Galactic-average, three-part power-law IMF i,(m) ex 
TO~"', where dn = ^(m) dm — ^^lI^^) dlogif^m is the num- 
ber of objects in the mass interval m,m + dm and logiow, 
logiom -I- dlogioTO, respectively. ■^lI'ti) = (rnlnlO) ^(to) is 
the "logarithmic IMF". The standard IMF has (Kroupa 
2001) 

ao = +0.3, 0.01 <m< O.OSMq; 



ai = +1.3, 0.08 <m< 0.50Mq; 
a2 = +2.3, 0.50 < m < 50Mq. 



(1) 



Initially all objects, stars and brown dwarfs, are assumed 
to be in binary systems with component masses chosen 
randomly from the IMF. The total binary proportion 



/tc 



Ni 



bin 



Nhin + Nsi 



where A'bin and A"sing are the number of binary and single- 
star systems respectively, is then equal to unity. The mod- 
els initially consist of 10"* stars and brown dwarfs in 5000 
binaries. 

Two models are constructed, both with a spherical 
Plummer density profile for the stars and the gas initially. 
The half-mass radius is -Ro.5 = 0.45 pc for model A with an 
initial central number density pc — lO**'* stars/pc'^, and 
i?o 5 = 0.21 pc for model B with pc = lO^'* stars/pc^. 
Both models are embedded in a gas potential with twice 
the mass in stars and with the same density profile as the 
stellar component. The gas is removed due to the action 
of the O stars after 0.6 Myr and the dynamical evolution 
is integrated for 150 Myr. By assuming a star-formation 
efficiency of 33 per cent and removing the residual gas 
within a thermal time-scale, the initially embedded clus- 
ter (the proto-ONC) expands rapidly, matches the ONC 
by about 1 Myr, and by 100 Myr it forms a bound star 
cluster which resembles the observed Pleiades and which 
forms the nucleus of an expanding association. 
The MF measured in the Pleiades today may bear witness 
of these events, because if the embedded cluster was signif- 
icantly mass segregated prior to gas removal, then mostly 
the low-mass stars will have been lost, and the Pleiades 
would then be deficient in very-low-mass stars. 
The evolution of the two cluster models is discussed in 
KAH and compared to observational constraints on the 
structure, kinematics and binary-star properties for both 
the ONC and the Pleiades cluster. The core radii as well 
as the number of stars within these radii predicted at 0.9 



Myr and lUU Myr by both models are m reasonable agree- 
ment with observed values for the ONC and the Pleiades 
respectively. Available velocity dispersion measurements 
are also consistent with the models. The observed ONC 
radial density profile is well reproduced by model B but 
not model A, whereas for the Pleiades the best fit is given 
by model A. Concerning the period distribution for late- 
type binaries, both models are consistent with the Pleiades 
data, although an intermediate model in term of central 
density would be better. For the ONC, only model B leads 
to agreement with the measured proportion of binaries 
with periods P > 10^ days, model A retaining too many 
binaries. Overall though, both models fit the ONC and 
the Pleiades data quite well, suggesting that the Pleiades 
cluster most probably formed from an ONC-like object. 
Therefore the initial state of the Pleiades will have been 
much more concentrated than presently observed with a 
half-mass radius between 0.2 and 0.5 pc. 
In this paper we extend this study by a comparison to 
the Pleiades MF which is now constrained down to the 
sub-stellar regime. The system MF which is constructed 
by counting only the system masses is compared with the 
observed MF, because the observational surveys cannot 
resolve binary stars. It is interesting to investigate the 
Pleiades MF to see if it may have changed and whether 
it differed from the standard IMF, given that the search 
for IMF variations constitutes one of the most important 
issues of star formation work. 



3. The observed Pleiades mass function 

Moraux et al. (2003) performed a deep and extended sur- 
vey of the Pleiades cluster using the CFII-12K camera. 
The survey covered 6.4 sq.deg. reaching up to 3 degrees 
from the cluster's centre and the covered mass range ex- 
tends from 0.03 to O.45M0. The authors used the NextGen 
(Baraffe et al. 1998) and Dusty (Chabrier et al. 2000) mod- 
els to convert their luminosity function into a mass func- 
tion and they used the Stauffer & Prosser open cluster 
database^ to complete their data for m > 0.4Mq. They 
thus computed the mass function from 0.030Mq to IOMq 
and found that it is reasonably well-fitted by a log-normal 
function in logarithmic miits (Moraux et al. 2003, Fig. 9) 
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with mo ~ O.25Af0 and <t — 0.52. It is shown on Fig.Q] 
However, the lowest mass point of the MF looks a bit 
different and is located much above the log-normal fit. 
Similar features are observed for other clusters around the 
same spectral type (M7-M8) and Dobbie et al. (2002) ar- 
gued this reflects a sharp local drop in the luminosity-mass 
(L-M) relationship, due to the onset of dust formation in 
the atmosphere around T^ff ~ 2500 K. A change in the 
slope of the L-M relationship implies that the masses of 
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objects with spectral types later than M7 may be sig- 
nificantly underestimated by the current NextGen and 
Dusty models and that the number of objects in this low- 
est mass bin may be overestimated. By applying the em- 
pirical magnitude-mass relation given by Dobbie et al. 
(2002; their Fig. 3) for the Pleiades low mass stars and 
brown dwarfs, we find that the lowest mass point of the 
Pleiades MF goes down whereas higher mass points go up 
(see Fig. nj. The log-normal fit does not change however 
and becomes even better. 
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Fig. 1. The Pleiades mass function. The open triangles 
are from the Prosser and Stauffer Open Cluster Database 
and the filled circles are the data points from Moraux et al. 
(2003). The open squares have been obtained by applying 
the empirical magnitude-mass relation given by Dobbie et 
al. (2002). The log- normal fit (eq.|2) found by Moraux et 
al. (2003) is shown as a solid fine. 



Even if this dust effect has still to be investigated, it sug- 
gests that the lowest mass point of the Pleiades MF may 
come down in the future and we consider the log-normal 
distribution as a good approximation of the observed MF 
in this mass domain. 

For m ^ 2Mq however, this fit does not seem to be steep 
enough. The 2-AMq MF points are more than 2-sigma 
below the log-normal fit. 



4. Comparison models/observations 

4.1. The mass function 

The system MFs at i = 100 Myr resulting from the dy- 
namical and stellar evolution of KAH's models, assuming 
random pairing and a binary fraction /tot = 1 at i = are 
compared to the log-normal fit of the observed Pleiades 
MF in Fig. 13 over the entire mass range. The system MF 
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Fig. 2. Comparison of the log- normal fit of the observed 
Pleiades MF (shown as a solid line) to the model MFs. The 
dashed (dotted) histogram represents the system mass 
function at an age of 100 Myr from model A (B) of KAH 
for all systems within R<15 pc of the model cluster center. 
They have been scaled to match the Pleiades data around 

TO ~ IMq. 

Both models reproduce well the log-normal distribution 
for TO ^0. IMq which indicates that the observed MF 
(data points in Fig. ^ is also steeper than the model MF 
for TO ^ 2 Mq. This discrepancy cannot be the effect of the 
stellar evolution since the models include its treatment but 
may indicate that the true stellar IMF in the Pleiades may 
have been steeper than Salpeter {a = 2.35) above 2Mq. 
It is useful to note that the KAH models do not constitute 
a fit but are a prediction of the Pleiades MF in the stellar 
and the sub-stellar regime under the assumption that the 
Pleiades IMF was identical to the standard IMF (eq. Q 
Kroupa 2001). Our result here suggests that this may in- 
deed have been the case except for the possibly steeper 
IMF for TO>2M0. 

For to^O.IMq down to the sub-stellar regime, model A 
yields a better fit than model B which suggests that the 
Pleiades may not have formed from an embedded cluster 
as concentrated as model B. The larger number of isolated 
BD systems predicted by model B compared to model A 
in Fig. 121 results from the early release of brown dwarf 
companions from their primary in the initially denser en- 
vironment of model B. The observed period-distribution 
of late-type Pleiades binaries seems also to be better re- 
produced with model A for P < 10^ days (see KAH, their 
Fig. 11), although an intermediate model in terms of cen- 
tral density would lead to improved agreement for longer 



periods, ihererore, our results tend to lavour modei A 
but model B cannot be excluded since it is still consistent 
with the observational result within the uncertainties (es- 
pecially for the lower mass bins where the effect of dust 
formation still needs to be investigated). 
By 100 Myr the binary fraction is not equal to one any 
longer and depends on mass. During the early stage of the 
cluster's life, many binaries are disrupted because of grav- 
itational interactions. The binding energy being weak for 
low mass binaries, many brown dwarfs (BDs) have been 
freed and the BD-BD binary fraction predicted by the 
models is ^ 20% (see KAH, their fig. 9) which is consistent 
with current observational results ranging from 15-50%. 
Martin et al. (2003) find - 15%±1^ for resolved Pleiades 
systems with separations from 7-12 AU but Pinfield et 
al. (2003) suggests that the Pleiades total (i.e. unresolved 
in general) binary fraction is from 30-50%. Overall, the 
observational constraints on late-type binaries are in rea- 
sonable agreement with the models, as shown by KAH. 
The binary proportion of the massive stars decreases first 
through disruption and then increases as new companions 
are captured. 

The mass function of all the objects counted individually 
is represented by a histogram in Fig. O for both models. 
It corresponds to the Pleiades MF we would observe if we 
were able to resolve all the binaries. If we compare this 
single star MF with the standard IMF indicated by the 
dots, we find that the overall shapes are similar, although 
model B shows a marginal deficit of BDs compared to 
model A. This indicates that the main difference between 
the Pleiades MF observed at an age of r^ 100 Myr and 
the Pleiades IMF is mainly due to the binarity effect and 
that the cluster dynamical evolution had little effect on 
its mass distribution. 

The Galactic disk log-normal MF obtained by Chabrier 
(2003, eq. 17) by fitting L and T dwarf samples (Chabrier 
2002) is found to be slighty above the histograms in the 
sub-stellar domain. This suggests that either the Pleiades 
IMF is different from the Chabrier's field IMF in this mass 
range or that the number of brown dwarfs in the field has 
been overestimated. The ratio 



3.5 



i? = 



NobjiO.02-0.08MQ) 
Nobj{0.l5-1MQ) 



where Nobj is the number of objects in the respective mass 
range is found to be equal to 1.07 ± 0.17 for Chabrier's 
Galactic disk log-normal MF, whereas it is equal to 0.81 ± 
0.05 for the three-part power-law standard Pleiades IMF. 

4.2. Mass segregation 

In Fig. 01 the fraction of systems per mass bin contained in 
the cluster centre (R<2 pc) and predicted by the models at 
t = 100 Myr is compared with observational constraints. 
The solid line corresponds to model A and the dashed line 
to the initially more concentrated model B. The observa- 
tional results for the Pleiades are indicated by the filled 
circles. The data points come from Pinfield et al. (1998) 
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Fig. 3. The single object mass function for model A 
(dashed histogram) and model B (dotted histogram) at 
t = 100 Myr and R<15 pc. The shaded region represents 
the field mass function from Chabrier (2003) with uncer- 
tainties. The dots correspond to the shape of the standard 
three-part power-law IMF. All models are scaled to agree 
between 0.5Mq and IMq. 



for the stellar domain and from Moraux et al. (2003) for 
the sub-stellar regime. The King density profile found by 
these authors for each mass bin has been de-projected, 
thus providing the spatial density which has then been 
integrated between and 2 pc to obtain the number of 
systems in the inner part of the cluster. 
Model predictions are consistent with observations. The 
more massive stars are concentrated around the cluster 
centre whereas the low mass objects are much more spread 
out throughout the cluster. 

That the mass segregation at i = 100 Myr is similar for the 
two models indicates the cluster to be completely relaxed 
at this age, and that the memory of the initial concen- 
tration has been lost. However, if we compare the mod- 
els after 1 Myr, mass segregation has not occurred yet 
for model A whereas it is already present for model B 
(see KAH, fig. 15). This model evolved more rapidly than 
model A by virtue of its larger initial density and thus 
shorter relaxation time. 

4.3. Evolution of the mass function 

In the previous sections it has been shown that the 
Pleiades cluster is dynamically evolved and that mass seg- 
regation has occurred, as is also found to be the case by 
Portegies Zwart et al. (2001). Here we consider the de- 
tailed structure of the MF to infer if the preferential loss 
of stars of a certain mass may be evident in the models. 
The ratio of the present-day MF at i = 100 Myr to the 
IMF is plotted for models A and B in Fig. [3 The compar- 
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being prcicrcntialiy lost, ihe enect is not very signihcant 
though, amounting to about 16 per cent. 
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Fig. 4. Ratio of the number of systems within R<2 pc to 
the total number of cluster systems (R<15 pc) per mass 
bin at t — 100 Myr. The sohd line corresponds to model 
A and the dashed line to model B. The stellar data points 
are from the observational survey of Pinfield et al. (1998) 
and the sub-stellar one is from Moraux et al.(2003) (filled 
circles). 

ison between the cluster population at 100 Myr and the 
initial cluster population shows that the evolved popula- 
tion is depleted above about AMq due to stellar evolu- 
tion. However there are no significant evolutionary effects 
that deplete the population of stars having a mass from 
2 — 4:Mq. Therefore the difference between the observed 
Pleiades mass function and the models in this mass range 
is more likely due to binarity or a different IMF. 
For model A the ratio is quite constant for less massive 
stars (and BDs) thus demonstrating that the same frac- 
tion of objects has been lost from the cluster indepen- 
dently of mass and thus that energy equipartition did not 
play a major role in de-populating the embedded clusters 
but rather the violent process given by rapid gas expulsion 
during which the model clusters loose about 2/3 of their 
population. 

The larger ratio for model B comes about because this 
model is initially more concentrated thus ultimately lead- 
ing to a more massive bound cluster because a larger frac- 
tion of the population evolves to a more bound state due 
to energy equipartition before gas expulsion, as already 
noted by KAH and considered in detail by Boily & Kroupa 
(2003a, 2003b). For model B the ratio decreases from 0.45 
at IM0 to 0.38 below 0.08 M©. This comes about in this 
model, which has a short equipartition time-scale, because 
less-massive objects end up with a lower cluster binding 
energy due to energy equipartition during the embedded 
phase. Gas expulsion then leads to the low-mass objects 
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Fig. 5. Fraction of individual objects remaining inside the 
cluster (R<15 pc) at i = 100 Myr compared to the initial 
number of objects at i = 0. The solid line corresponds to 
model A and the dashed line to model B. The ratio of the 
individual object MFs is shown, not of the system MFs. 



5. Conclusions 

The Pleiades is a young (« 120 Myr) close-by (« 120 pc) 
Galactic cluster and therefore allows very detailed obser- 
vational scrutiny. The available data together with the- 
oretical work suggest that the Pleiades may have been 
born similar to the ONC, and that its subsequent evolu- 
tion passed a violent phase as a result of rapid gas expul- 
sion once the massive cluster stars provided sufficient feed- 
back to the young nebula. The theoretical models suggest 
that initially highly concentrated states (model B), with 
Pc ~ 10^-* stars/pc^, evolve to significant mass segrega- 
tion prior to gas expulsion, while less-concentrated models 
with Pc ~ lO**-^ stars/pc^ (model A) do not. As a result, 
the former looses slightly more of its low-mass members 
relative to higher-mass ones than the latter (Fig. El . 
Both models are in good agreement with the observed MF 
in the Pleiades (Fig. |21l, although model A yields a better 
fit which may suggest that the Pleiades precursor had a 
concentration near pc ~ lO**'* stars/pc^ rather than being 
more concentrated. The cluster's stellar IMF could thus 
have been very similar to the Galactic-field average IMF 
for stars with m^2 Mq. No significant stellar IMF varia- 
tion can thus be detected, in comparison to the Galactic- 
field. The sub-stellar part of the Pleiades mass function 



also seems correctly reproduced by the dynamical evolu- 
tion of dense clusters. This suggests only marginal evap- 
oration of the lowest cluster members relative to higher 
mass stars on a timescale of 100 Myr. 
The Pleiades do appear to have a deficit of stars with 
m^2 Mq (Fig. [21 compared to the number expected from 
the average IMF indicating that the Pleiades IMF may 
have been steeper, a^ > 2.3 for m > 2Mq. This is inter- 
esting, because there exists evidence that the stellar IMF 
corrected for binarity for massive stars is probably steeper 
than a Salpeter IMF (Sagar & Richtler 1991; Kroupa & 
Weidner 2003). These authors found a^ « 2.7 whereas an 
apparent Salpeter IMF indicates aapp — 2.35. However, 
it will be necessary to compute additional models with 
different values for the star-formation efficiency and mass 
to map-out feasible Pleiades precursors and the range of 
physical conditions within the solution space. For exam- 
ple, a lower initial density and thus a presumably larger 
loss of massive stars during gas expulsion may reduce the 
discrepancy found here of the measured MF and the stan- 
dard IMF. 

The Pleiades show very pronounced mass segregation 
(Fig. 21) which is well reproduced by both models. This in- 
dicates that the mass segregation that develops in model B 
prior to gas expulsion is forgotten by 100 Myr. The 
Pleiades is therefore in a completely relaxed state, hav- 
ing almost completely lost memory of its initial condi- 
tion. These, however, are reflected in the properties of its 
binary population, because the initial concentration lim- 
its the binary-star periods that survive to old cluster age 
(KAH). 
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